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■ Abstract. The effects of the Galactic bar on the velocity distribution 

fNJ , of old disc stars in the Solar neighbourhood are investigated using high- 

resolution 2D test particle simulations. A detailed orbital analysis reveals 
that the structure of the U — V distribution in these simulations is closely 
related to the phase-space extent of regular and chaotic orbits. At low 
angular momentum and for a sufficiently strong bar, stars mainly follow 
chaotic orbits which may cross the corotation radius and the U — V con- 
tours follow lines of constant Jacobi's integral except near the regions 
occupied by weakly populated eccentric regular orbits. These properties 
can naturally account for the observed outward motion of the local high 
^ | asymmetric drift stars. 

o 
o 
o 

| 1. Introduction 

s: 

Q-f The observed local stellar velocity distribution in the U — V plane, where U is 

the radial velocity positive towards the anti-centre and V the azimuthal velocity 
positive towards galactic rotation, betrays a large stream of old disc stars with 
an asymmetric drift S ~ 40 — 50 kms -1 and U > (e.g. Dehnen 1998), which 
hereafter will be referred to as the "Herculis" stream according to a co-moving 
Eggen group. One of the most reasonable interpretation of this stream relies on 
the influence of the Galactic bar, which is believed to have a corotation radius 
Rcr ~ 3.5 — 5 kpc and a major-axis presently leading by 99 as 15° — 45° the 
Sun's azimuth. Raboud et al. (1998) have shown that iV-body models of the 
Milky Way do produce a mean outward motion of disc particles at realistic posi- 
tions for the Sun relative to the bar and suggest that the Herculis stream involves 
stars on hot orbits, i.e. orbits with a high enough value of Jacobi's integral to 
cross the corotation radius, whereas Dehnen (2000, hereafter D2000) relates this 
velocity mode and the main low-asymmetric drift mode of the observed U — V 
distribution to the coexistence near the outer Lindblad resonance (OLR) of the 
inner anti-bar and outer bar elongated periodic orbits replacing the circular orbit 
close to the OLR in a rotating barred potential (i.e. the same idea introduced 
by Kalnajs 1991 to explain the Hyades and Sirius streams), and attributes the 
valley at S ~ 30 kms -1 between the two modes (respectively "OLR" and "LSR" 
mode in his terminology) to stars on unstable OLR orbits. This paper analyses 
the role of bar induced chaotic orbits in this context and confronts both alterna- 
tives using a similar 2D modelling technique as in Dehnen's work. The details 
are presented in Fux (2000). 
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2. Resonances and periodic orbits 

The working potential is the same as in D2000, consisting of the sum of an 
axisymmetric logarithmic potential with constant circular velocity V and of a 
rotating quadrupole term with a bar strength F = 0.15, where F is defined 
as the maximum azimuthal force on the circle passing through the ends of the 
bar normalised by the axisymmetric radial force at same radius (F = 8.89a in 
Dehnen's notation). The position in space will be parametrised by the galac- 
tocentric distance relative to the OLR radius -Rolr = (1 + 1/v / 2)-Rcr an d the 
angle ip relative to the major axis of bar. A distance of i?/i?oLR = 1-1 implies 
R C r = 4.26 kpc if R = R Q = 8 kpc. 

Figure 1 (left) shows the characteristic curves in the Hamiltonian— x plane 
(at y = and y > 0) of the main periodic orbits existing near the OLR and 
closing after one rotation in the rotating frame, using the labelling of Contopou- 
los & Grosbol (1989). The curves essentially resemble those of the circular orbit 
and of the resonant orbits in the axisymmetric limit, except near the resonances 
of the circular orbit where bifurcations may transform into gaps, as is the case 
at the OLR (e.g. Sellwood & Wilkinson 1993). Figure 1 also displays some 
members of the xi(l) and xi(2) families. The thick orbits are examples of the 
perpendicularly oriented orbits that replace the circular orbit near the OLR and 
which are responsible for the bimodality in the observed U — V distribution ac- 
cording to D2000. In both families, the orbits develop loops at high eccentricity, 
i.e. on the part of the characteristic curves emerging leftwards from the zero ve- 
locity curve and where x ^ 0, and the x±(l) orbits are stable over a remarkable 
eccentricity range. Figure 2 shows the location in the U — V plane of the —4/1, 
—2/1 (OLR) and —1/1 resonances in the axisymmetric part of the potential 
for different space positions, and of the periodic orbits passing through these 
positions in the total axisymmetric+bar potential. At <p = 90° most resonant 
periodic orbits are stable, except for instance the x\{2) orbit at i?/i?oLR ^ 1-1, 
whereas at cp = 0, most of them are unstable. 

3. Hot orbits 

In a rotating barred potential, the value of the Hamiltonian H = ^x + <b e g(x), 
where $ e g(x) is the effective potential, is the only classical integral of motion, 
known as Jacobi's integral. Two critical values of this integral are related to 
the Lagrangian points L1/2 and £4/5, i.e. respectively the saddle points and 
maxima of 3> e ff(x). If H > H\2 = <& e s(Li/2), an orbit is susceptible to cross the 
corotation radius, defining the hot orbit category, and if H > H45 = <& c fj(L 4 / 5 ), 
an orbit is no longer spatially confined by this integral. 

If U and V are measured relative to the galactic centre, the contours of 
constant H in the U — V plane are circles centred on (V, U) = (Rftp, 0), where 
Qp is the pattern speed of the bar, and of radius \/2{H — <& e ff)- Figure 2 gives the 
H\2 and contours as a function of space position for the working potential. 
The hot orbits are those outside (low V) the former contour. At R = i?oLR and 
ip = 30°, the H\2 contour crosses the V-axis at S = V —V pa 0.17714 = 35 km s _1 
for V = 200 kms^ 1 , placing the Herculis stream in the hot orbit zone. 
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Figure 1. Left: characteristic diagram of the relevant periodic orbit 
families for F = 0.15. The dashed line indicates the zero velocity curve 
and -ffoLR is the Hamiltonian of the circular orbit at the OLR in the 
axisymmetric part of the potential. Middle and right: sequences of 
xi(l) and xi(2) orbits, with the thick line representing the orbit of 
lowest apocentre and of lowest Hamiltonian respectively. The orienta- 
tion of the bar is sketched by the shaded ellipse. In all frames, full and 
dotted lines stand for stable and unstable orbits respectively. 



4. Regular and chaotic regions of phase-space 

The next step after the periodic and hot orbits is to determine the extension 
in phase-space of the regular orbits trapped around the stable periodic orbits 
and of the chaotic regions. A convenient tool for this purpose are the Liapunov 
exponents, which describe the mean exponential rate of divergence of two trajec- 
tories with infinitesimally different initial conditions. The largest such exponent, 
Ai, determines the degree of stochasticity of the orbits and vanishes for regular 
orbits. Figure 2 gives the divergence timescale t\ = 1/Ai of the orbits as a 
function of planar velocity and space position. 

The Liapunov divergence timescale can be as low as a few orbital times for 
chaotic orbits. As expected, the stable periodic orbits always lie within reg- 
ular regions and the unstable orbits within chaotic regions. Furthermore, the 
H\2 contour appears as the average transition from regular motion at low H to 
chaotic motion at high H, and thus most hot orbits are chaotic. Due to the 
four-fold symmetry of the potential, the diagrams at ip = and ip = 90° are 
symmetric with respect to U = 0, and more generally, diagrams at supplemen- 
tary angles are anti-symmetric to each other in U. At ip = 90°, the —2/1 and 
—1/1 resonance curves are mainly embedded within broad arcs of regular orbits 
spaced by chaotic regions at large H, while at <p = 0, these regular arcs occur 
between the resonance curves. At intermediate angles, the regular and chaotic 
regions become offset from the resonance curves and the [/-symmetry brakes. 
In particular, for <p ~ 30° and for the displayed radial range, a prominent region 
of regular orbits trapped around stable eccentric xi(l) orbits extends down to 
negative U, bounded roughly by the OLR curve on the right and penetrating 
well inside the hot orbit zone, whereas the positive U part of the OLR curve is 
surrounded by a wide chaotic region extending somewhat inside the H\i con- 
tour and coinciding very well with the observed U — V location of the Herculis 
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Figure 2. Liapunov divergence timescale of the orbits in the U — V 
plane as a function of space position, for F = 0.15. The horizontal and 
vertical axes of each frame are V/V and U/V respectively, with the 
origin at the circular orbit of the underlying axisymmetric potential. 
The timescale is in units of local circular period tn- The dark and white 
regions respectively indicate regular and chaotic orbits, the various full 
and empty symbols the stable and unstable periodic orbits, the circular 
arcs open towards the right the H\i and contours, and the dash- 
dotted, solid and dashed curves the —1/1, —2/1 and —4/1 resonances 
in the axisymmetric part of the potential. 

stream. Note that since 2D axisymmetric potentials are always integrable, all 
chaotic regions reported here owe uniquely to the bar. 

A detailed investigation reveals that the regular orbits associated with the 
stable low eccentricity periodic orbits of the xi(2) family exist for R/Rq,lr ^ 11 
and over an angle range around <p = 90° increasing as R decreases. Elsewhere, 
unstable x\(2) orbits have completely substituted to the regular x\{2) islands. 
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Figure 3. Velocity distribution as a function of bar parameters in the 
test particle simulations with F = 0.15. The distributions are averaged 
over the time interval going from 55 to 65 bar rotations. The axes of 
the frames and the resonance and H contours are as in Figure 2. 

In particular, there is no x\{2) streaming motion expected at R/Rolr = 1-1 
and ip = 30° in the present working potential. 

5. Phase-space crowding 

Finally, we want to know how nature populates the available orbits. This is done 
by test particle simulations starting with exactly the same initial axisymmetric 
distribution functions f and growing the bar the same way (i.e. over 2 bar 
rotations) as in Dehnen's default simulations. However, we do not apply the 
claimed high-resolution and noise-free backward integration technique based on 
the conservation of the phase-space density along orbits. The reason is because 
much longer integration times than the few bar rotations adopted in D2000 are 
necessary to get phase-mixed and quasi-steady evolved distribution functions, 
and that for integration times larger than about 10 bar rotations, the fine-grained 
U — V distributions resulting from this technique become highly fluctuating on 
small scales and even noisy in the chaotic regions and thus must be smoothed 
to yield the physical coarse-grained distributions. Therefore, we decided to use 
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the standard forward integration technique, sampling f a with N = 10 6 particles, 
and to average the velocity distributions, derived within finite volumes of radius 
R/80, over 10 bar rotations to improve the statistics and also to minimise the 
phase-mixing problem. The results are depicted in Figure 3. Here the distance 
parameters in f Q always scale as R and the truly changing parameters are those 
of the bar. Hence frames at different R/Rolr come from distinct simulations. 

The U — V distributions seem to present the same symmetry properties in U 
as described in the previous section for the Liapunov divergence timescale. This 
is obviously not the case in D2000, as a consequence of the too short integration. 
The high angular momentum peak of the velocity distributions also coincides 
better with the trace of the non-resonant x±(l) orbits at U = 0. The regular 
low eccentricity x±(2) regions also produce a peak in the velocity distributions, 
which always lies inside the H\2 contour. 

In the hot orbit zone, the chaotic regions appear to be more heavily crowded 
than the regular regions. This is because the phase-space density in regular re- 
gions remains roughly constant during the simulations, whereas there is a net 
migration of particles on hot chaotic orbits from the inner to the outer space 
regions, and because the chaotic and regular regions are decoupled from each 
other, i.e. chaotic orbits cannot visit the regions occupied by regular orbits 
and vice versa. Hence at large H, the velocity distributions are dominated by 
chaotic orbits which are forced to avoid less populated regular regions. In the pa- 
rameter range considered here, these forbidden regions are essentially the broad 
regular orbit regions around the stable eccentric xi(l) orbits. In particular, at 
R/Rolr = 1-1 and if = 30°, this x\(l) region lies at U <; 0, as mentioned 
above, and thus the main overdensity of hot chaotic orbits occurs at positive U, 
providing a natural interpretation of the Herculis stream. The depression on the 
right of this stream may be attributed to the decline of the hot orbit population 
as H — ► H\2- The high eccentricity xi(l) region is part of Dehnen's valley, but 
represents regular OLR orbits, and the unstable x*(2) orbit does not always fall 
within this valley (e.g. at R/Rolr = 1-2 and (p = 30°). 

Another property of the hot chaotic regions is that the contours of the veloc- 
ity distribution follow the lines of constant Hamiltonian. This can be understood 
in terms of the Jeans theorem (see also the answer to Dehnen's question): in 
two-dimensions, chaotic orbits have no other integral of motion than Jacobi's 
integral and thus the distribution function in their neighbourhood depends only 
on this integral, i.e. on the value of the Hamiltonian. 
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Questions 

Ortwin Gerhard: The spiral arms outside the Milky Way's bar region must 
modify the population seen in the U — V diagram. Could you say how? 

Roger Fux: The effect of spiral arms is not yet very clear. As evolving transient 
features, they certainly introduce a time dependency of the U — V distributions 
and probably increase the amount of chaos. Since Jacobi's integral is then no 
longer conserved, they may also introduce some diffusion from the regular regions 
to the chaotic regions and vice versa (and the same remark is true for massive 
compact objects like giant molecular clouds). iV-body simulations displaying 
spiral arms reveal time averaged U—V distributions that are qualitatively similar 
to those of the test particle simulations presented here, but with an azimuthal 
phase shift apparently related to the twist of the potential well induced by the 
spiral arms outside the bar. 

Walter Dehnen: How can stars stay on chaotic orbits, which due to their H 
are not bound to the galaxy? 

Roger Fux: The timescale for chaotic hot orbits with moderate H to escape 
the galaxy is much larger than the orbital timescale and even the Hubble time 
(see for instance Kaufmann & Contopoulos 1996, A&A 309, 381). In particular, 
once phase-mixed, the distribution function in the neighbourhood of such orbits 
may be considered as almost steady, which is a necessary condition to apply the 
Jeans theorem. 



